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Abstract 

The production in high-energy hadron collisions of a pair of jets with large rapid¬ 
ity separation is studied in an improved BFKL formalism. By recasting the analytic 
solution of the BFKL equation as an explicit order-by-order sum over emitted gluons, 
the effects of phase space constraints and the running coupling are studied. Particular 
attention is paid to the azimuthal angle decorrelation of the jet pair. The inclusion of 
sub-leading effects significantly improves the agreement between the theoretical pre¬ 
dictions and recent preliminary measurements from the D0 collaboration. 


1 Introduction 


Fixed-order, renormalization-group-improved QCD perturbation theory has been remarkably 
successful in describing jet physics at high-energy colliders. For example, the bulk of the 
inclusive jet distribution at the Tevatron pp collider is well described by folding subprocess 
cross sections calculated at NLO with parton distributions extracted from deep inelastic 
scattering. However there are certain situations where such an approach may be expected to 
fail. If in the production of a pair of large Et jets the separation in rapidity /S.y of the jets 
becomes large, then higher-order perturbative corrections become more and more important. 
Essentially, for Ay S> 1 the subprocess cross section has an expansion in powers of agAy 
rather than 

Dijet production with a large rapidity separation is an example of a ‘two-large-scale’ pro¬ 
cess in perturbative QCD, where large logarithms, in this case ln(s/E|,) ~ Ay ^ 1, arising 
from real and virtual soft gluon emission compensate the strong coupling in the perturbation 
series. These logarithms can be resummed using the techniques of Balitsky, Fadin, Kuraev 
and Lipatov (BFKL) |]^. Indeed it was first pointed out by Mueller and Navelet that dijet 
production in high-energy hadron-hadron collisions would be a particularly clean environ¬ 
ment in which to look for evidence of such resummation. They showed that the subprocess 
cross section was expected to increase at asymptotic separations according to d ~ exp {XAy) 
with A = asl21n2/7r. These ideas were taken further in Refs. [|], §, ^ |^. In particular, it 
was shown that the azimuthal decorrelation of the jet pair, resulting from multiple soft gluon 
emission in the rapidity interval between them, provides a particularly distinctive signature 
of BFKL dynamics, free of potentially confounding effects from dependence on parton dis¬ 
tributions. The leading BFKL resummation provides a quantitative prediction for the rate 
of decorrelation with increasing separation Ay. 

The azimuthal decorrelation has recently been measured by the DO collaboration at 
the Tevatron pp collider. The results are intriguing. The observed rate of decorrelation is 
larger than that predicted by an ‘exact’ next-to-leading order calculation based on the 2 —>■ 2 
and 2 —3 matrix elements (as implemented in the JETRAD program 0), but smaller than 
that predicted by both the leading BFKL resummation p, ^ and an ‘improved’ BFKL 
calculation incorporating certain subleading kinematical effects 0. In fact the data agree 
best with the predictions of the HERWIG parton-shower Monte Carlo [Q, based on DGLAP 
multigluon emission with angular ordering. 

The bulk of the theoretical attention concerning BFKL dynamics has focused on its 
application to HERA ep collider physics - originally, to the behavior of the structure function 


F 2 at small x, and latterly to more exclusive quantities (for a review see Refs. [0,0 
references therein). In particular, ‘forward jet’ production has been studied in Refs. 
131, [0, ^ |T6|, and forward single particle and single photon production in Ref. 


and 


12 


The 

possibility of observing ‘BFKL gluons’ in small-a: deep inelastic scattering has been discussed 
in Ref. [j0|. One result of this activity is the realization that subleading corrections to the 
leading BFKL resummation, from effects such as the running coupling and phase space, are 
likely to be numerically important in practice; see for example Ref. [|0|. Unfortunately the 
complete perturbative next-to-leading logarithmic corrections are not yet available in a form 
which allows them to be incorporated in a phenomenological analysis, although there has 
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been considerable theoretical progress towards this goal 

The aim of the present work is to study such subleading effects and to see whether they 
could account for the discrepancy between the BFKL predictions and the Tevatron data. We 
do this by recasting the original analytic resummed expression for the parton cross section [Q 
in the form of an event generator. This not only allows kinematic subasymptotic constraints 
(energy conservation, for example) and other corrections to be readily implemented, but 
also gives a better overall picture of the events with large dijet rapidity separation (for 
example, the distribution of the accompanying BFKL ‘minijets’). A further advantage is 
that experimental acceptance cuts on the jets are straightforward to impose. A similar 
Monte Carlo approach was used recently in the fixed coupling limit to look at transverse 
energy flow in dijet production. 

In this paper we concentrate on several issues. First, we set up the calculational frame¬ 
work and demonstrate that we are able to reproduce the results based on the analytic 
resummed expressions of Refs. il. 1- Second, we examine the effect on the predicted 
azimuthal decorrelation of including various subasymptotic effects, including kinematic con¬ 
straints and the running coupling constant. We will show that these effects are indeed 
important in the kinematic domain accessible to the Tevatron experiments. More general 
applications of our formalism, together with a more detailed comparison with data, will be 
presented elsewhere. 


2 BFKL formalism for dijet production 

2.1 Summary of analytic results 

We wish to describe events in hadron collisions containing two jets with relatively small 
transverse momenta PtitPt 2 and large rapidity separation Ay = yi — y 2 . In the limit where 
the minimum jet transverse momentum Pt is small compared to the parton-parton center- 
of-mass energy and Ay is large, the inclusive dijet cross section can be written analytically 
to leading logarithmic order in the BFKL approach [Q. 

For jets produced in gluon-gluon collisions (the qq initial state can be included via the 
effective subprocess approximation; see below) the differential cross section is given by 


da, 
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.2 ^^2 0^3 ^3 9 ^ ^ - 2 tT . 
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dpl^dpl^dAcj) 2p^y^2 
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dz exp (2tXn{z) + izln(p|i/p|2)) (1) 


with t = asCAAy/iT and 


Xn{z) = Re ['0(1) - "0 (^1 + 1^1) + 


Here ip is the logarithmic derivative of the gamma function. We have dehned 

A0 = 101 - 02| - vr 

so that A0 = 0 when the two jets are back-to-back in the transverse plane. 


( 2 ) 

( 3 ) 
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Integrating over the dijet transverse momenta above a fixed threshold Pt then gives for 
the azimuthal distribution 


da, 
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Pti ^Pt2^^T 
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The total subprocess cross section is simply proportional to Co(t): 




TT 
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Its asymptotic behavior is determined by 

= 1 
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for f = 0 
for t ^ oo 


(4) 
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from which we see the characteristic BFKL prediction of an exponential increase in the cross 
section with large Ay (equivalently, large t). It can also be seen from that the average 
cosine of the azimuthal angle difference A0 defined above is proportional to Ci{t). In fact 
we have 

and as we shall see below, this falls off with increasing t, exhibiting the decorrelation expected 
with the emission of gluons in the rapidity interval between the jets. 

In what follows we will compare the analytic results just described with results from our 
formalism described in the next subsection. We will pay particular attention to the total 
subprocess cross section and the azimuthal angle correlation (||). 


2.2 Iterated solution for use in event generators 

The assumptions built into the leading order BFKL formalism that allow us to obtain an 
analytic solution for the cross section are not all easily satisfied in an experimental situation. 
In particular, implicit in the above solution are integrations over the transverse momenta of 
intermediate radiated gluons that extend to infinity. Furthermore nonleading effects (which 
lead for example to the running of as) are also neglected. And finally, the analytic solution 
(^,^ is symmetric with respect to the two observed jets, while experimental cuts may not 
be. 

In this section we solve the BFKL equation by iteration, which allows us to avoid these 
assumptions and obtain a solution that is more directly amenable to comparison with exper¬ 
iment. This solution amounts to ‘unfolding’ the summation over the intermediate radiated 
gluons and making their contributions explicit. It is then straightforward to implement this 
iterated solution in an event generator, as we describe below. 
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To obtain the iterated solution we begin with the differential cross section 
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d‘^PTid‘^PT2dAy PtiPt2 


f{PTl,PT 2 , Ay). 


( 9 ) 


The Laplace transform / of the function / with respect to Ay satishes the BFKL equation 
Dehning 

f{PTl,PT2,‘^)= / dA e f{pTl,PT 2 ,A) , 


( 10 ) 


where for convenience we use A = Ay , we have 

^f{PTuPT2,Uj) = d(pn -Pt2)^(01 - 02) + 

“ d^gT 
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TT^ 
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Qt 


f{PTl + qT,PT2,^) 


PTlfiPTl,PT2,Uj) 

Qt + (Pti + 


( 11 ) 


which is the BFKL equation for dijet production in hadron collisions. Infrared divergences 
from real gluon emission - the first term in the integral on the right-hand side - are cancelled 
by the virtual gluon contribution in the second term. Note that in writing Eq. (O) we have 
assumed that is fixed. 

Instead of solving (^) analytically, which would lead to the results shown above, we solve 
iteratively, using a slightly modihed form of the equation. Following Ref. |]^, we note that 
very low energy gluons are not resolvable and we therefore separate the real gluon integral 
into ‘resolved’ and ‘unresolved’ contributions, according to whether they he above or below 
a small transverse energy scale p. The scale p is assumed to be small compared to the other 
relevant scales in the problem (the minimum transverse momentum Pt of the ‘external’ jets, 
for example). We then combine the virtual and unresolved contributions into a single, hnite 
integral. The BFKL equation then becomes 


^f{PTl,PT2,Uj) = S{Pti -Pt2)^(01 - 02) + 


a,,C/ 


TT^ 


-^fiPTl + qT,PT2,^^) 
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( 12 ) 


The combined unresolved/virtual integral can be simplihed by noting that qx -C pri by 
construction because p^ -C and pxi > Pt- Therefore 


/(PTi + qr.Pn,^) ~ f(PTi,PT2,‘-j), 


(13) 


which allows us to write 

(a; -a;o)/(pTi,PT 2 ,a;) 


— 5{pxi Pt2)^(01 02) 


+ 
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fiPTl + qT,PT2,Uj), 
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where we have dehned 


UJq = 


/ (y,gC \ 
) 
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Pti 

4 + (Pti + QtY 


(15) 


The virtual and unr esolved contributions are now contained in ujq and we are left with an 
integral over resolved real gluons. 

We now solve Eq. (^) iteratively, and performing the inverse transform we have 


OO 

/(Pti,Pt2,A) = ^A) 


n=0 


(16) 


where the exact form of depends on whether the coupling ag is taken to be hxed or 
running. Strictly speaking, the derivation above only applies for hxed coupling because 
we have left ag outside the integrals. The modihcations necessary to account for running 
coupling are discussed below. 


2.2.1 Solution for fixed ag 


For the case where ag is hxed (evaluated at the scale P^, for example), the integral in ([^) 
is straightforward to evaluate and we obtain 


UJq — 


agC i 


TT^ 
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(17) 


Substituting into (p)^ , solving by iteration and then performing the inverse Laplace trans¬ 
form, we hnd 
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(19) 


The inverse transform has given rise to a set of nested integrals over the variables IZ, which 
can be interpreted as the rapidities of the emitted gluons. The diherential subprocess cross 
section is then given by 
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d‘^PTid?PT2d^y 
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The subprocess cross section is now expressed as an explicit sum over radiated gluons, with 
corresponding qxi and Yi integrals over their phase space. It is straightforward to implement 
this in a Monte Carlo event generator ,|^ and to impose energy conservation and experimental 
cuts via limits on the integration; see below for numerical results. 


An interesting feature of the results (1^) and (W) is the presence of the form factors 


ill 

Pti 


OLsCa^It^ 


and 


(PTI + YJi=l QTi)‘ 

(PTI + Er=l ^Ti)" 


XYn j'K 


respectively. These form factors arise from the resummation of the unresolved < [j? soft 
gluon emission in the rapidity interval A. In particular we see the modihcation of the naive 
zeroth-order perturbative result ^(pti +PT 2 ) by 


P 

Pti 


OsCaA/tT 


< 1 


for > fi"^, A > 0. 


( 21 ) 


This is a consequence of the fact that the emission of soft gluons reduces the probability of 
the dijets having equal transverse momenta and being back-to-back in azimuth. For A = 0 
all radiation is suppressed (in this approximation) and the form factor is equal to unity. 

We can make contact with the analytic results of the previous section by noting that the 
only additional approximation we have made is (|T^, which we used in the computation of 
(Uq, thereby neglecting momenta smaller than fi compared to pti- We therefore expect our 
result to agree with the analytic one up to corrections of /p^j). In fact one can show 
that the dijet cross section (PO]) integrated over the external transverse momenta 

= [ (fpTid‘^PT2 5(101 - 02| - vr - A0) 0(p^i - Pt)0{p12 - Pt) -T 2 - 1^2^^ WA 

corresponds exactly to the analytic result (|^) in the limit ^ 0. This will be illustrated 
numerically below. In the simulations based on Tevatron kinematics to be described in the 
following section, where pTi > 0(10 GeV), we will use values p = 0(1 GeV) such that the 
hnite /p^^) corrections are negligible. 


2.2.2 Solution for running 

Higher order corrections are known to lead to the running of the coupling constant 
0. They can therefore be taken into account by including momentum dependence in the 
coupling associated with the emission of each gluon. We shall see that this leads to two 
types of modihcation: (i) the factors of associated with resolved real gluon emission are 
simply evaluated at the scale of the emitted gluon; and (ii) the form factors associated with 
the unresolved real and virtual gluon emission get modihed slightly. 

We will include the running to lowest order, taking 


«s(g^) = 


/ A?) 


(23) 


^See also Ref. for a similar approach. 
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where 


b = 


33 - 2N- 
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127r 


(24) 


with Nf = 4. We then pull the factors of as in ( p!4D and ([T5| ) inside the integrals and make the 
substitution —> as{qT)- In addition, we must regulate the behavior of to prevent 

its becoming unphysically large as becomes small. This can happen for example in the 
momentum integration in the expression for ujq. Here we simply assume that the value of as 
freezes out below some scale Qq > A, i.e., we take 


as{q‘^ 


= ao = asiQl) < Ql 


= as{q ) 


for q^ > Qi 


(25) 


In practice choosing values A < Qq < fi so that ao is of 0(1), as we do below, gives results 
that are insensitive to the exact choice. 

The iterated solution to the BFKL equation in the running coupling case is then given 
by 
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Thus the result for running is obtained by the replacements 
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_1 
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in the hxed coupling results (|T^,^) for and and 

o^s(Pti)o^s(Pt2) 


a^ 


(29) 


(30) 


in the differential cross section (^). Note that as is equal to ag unless it is evaluated at 
a momentum smaller than the ‘freeze-out’ scale Qq. We maintain the hierarchy of scales 
A < Qq < /i <C Pt- 
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2.3 Cross section 


The calculation is completed by weighting the integrand in the dehnition of (Eqs. (P!U| , PBD ) 
in the subprocess cross section with parton distributions G(xi, Q^)G(x 2 , where, using 
the ‘effective subprocess approximation’, 

G{x,Q^) = g{x,Q^) + ^ + q{x,Q^)) . (31) 

q=u,d,s,c 

The parton momentum fractions Xi and X 2 are determined by the invariant mass and 
rapidity Y of the multijet hnal state: 


xi = y^e^ = -^e + Pt2 + 

X2 = (pTie~^ + pt2 + (32) 


where ^/s is the total collision energy. In the numerical studies to be described below we 
choose the factorization scale = P^. 

Note that the requirement xi,X 2 < 1 effectively imposes an upper limit on the trans¬ 
verse momentum {gri) integrals. This in turn means that the analytic results (^,^ are not 
reproduced in the presence of such a constraint, since they require the internal transverse 
momenta integrals to extend to infinity. In the original work of Ref. 0] (see also |^) the 
strong ordering of the rapidities was used to approximate the right-hand side of (|3^ by 


Pti 

Ts 

PT2 A/2-y 




(33) 


so that there was no longer any constraint on the qxi- An improved approximation where 
account was taken of the additional energy required for the multigluon emission was studied 
in Ref. |^. In our approach all kinematic constraints are applied directly to the multijet 
hnal state. 


3 Numerical Results 

In this section we present numerical results from our event generator for dijet production 
in the BFKL approach. We will take a brief look at the properties of the subprocess cross 
section a and then present predictions for the Tevatron pp collider with center-of-mass energy 
1.8 TeV. 

In our numerical computations we make the following choices for the relevant parameters. 
The minimum transverse momentum for each jet of the pair is Pt = 20 GeV, and for 
simplicity the two jets are assumed to have equal and opposite rapidities: pi = —y 2 = ^/2. 
Unless otherwise noted the scale p that dehnes the boundary between resolved and unresolved 


gluons we take to be /i = 1 GeV. In the fixed coupling case we evaluate as at the scale P|; 
for A = 200 MeV (as dictated by our choice of parton distributions [^]; see below) and our 
choice of Pt, this gives as = 0.164. In the running coupling case we again have A = 200 MeV 
and we choose Qq = 0.425 GeV so that Oq = 1. Finally, in our Tevatron calculations we use 
the leading-order parton distribution functions of Ref. 
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We begin with the subprocess cross section, shown in Figure 1(a) as a function of the 
dijet rapidity difference A, normalized to its value at A = 0. The results for fixed as are 
shown as open circles (the error bars shown are from Monte Garlo statistics), and we see 
the characteristic BFKL exponential rise with increasing A. We also see that the analytic 
result, shown as a solid curve, is well-reproduced by the iterated solution. The cross section 
for running ag is shown as the points marked ‘x’, and we see that the running of the coupling 
has the effect of slightly suppressing the increase in a with A compared to the fixed case. 
In Figure 1(b) we show (cos A0) computed from the subprocess cross section a for the same 
cases as in Fig. 1(a). The fixed-a* results (circles) agree very well with the analytic prediction 
(solid curve), while the running ag results show a slightly slower decorrelation. We will see 
a further slowing of the decorrelation when conservation of energy and parton distributions 
are included below. 

In Figure 2 we investigate the sensitivity of our results to the choice of gluon resolution 
cutoff fi. As noted above, we expect corrections to be of order so that the fi dependence 

becomes stronger as p approaches Pt, but for fi small enough there should be relatively little 
sensitivity. This is illustrated in Fig. 2 where we show the cross section a as a function of 
fi for two values of the dijet rapidity difference A. In this and the following figures we 
show results for running We see the expected behavior and also that larger values of 
A exhibit more sensitivity to fi. This is related to the fact, as we will see explicitly below, 
that larger A means more emitted gluons. This in turn leads to more sensitivity to the 
resolution parameter. We also note that the presence of other scales in the problem besides 
Pt (A, Qo) mean that arguments about the size of corrections are not rigorous. In any case, 
it is clear from the figure that our choice of /i = 1 GeV is sufficient to guarantee that we are 
not senstive to its exact value. 

We examine in Figure 3 how the cross section is distributed amongst the contributions 
from different numbers of emitted gluons for the same values of A as in the previous figure. 
As expected, larger rapidity differences between the jets allow for larger numbers of emitted 
gluons, so that the cross section peaks at higher Ug with increasing A. Note that the 
increasing area of the histograms simply reflects the rise of the cross section with increasing 
A. We should point out that there is not a. direct correspondence between emitted gluons and 
physical jets, so that these results for Ug should not be interpreted directly as a prediction 
for numbers of jets. By the same token, we note that the distribution in Ug depends to some 
extent on the resolution cutoff parameter fi — larger values of fi lead to fewer resolvable 
gluons. 

We now turn from the subprocess cross section to the total cross section integrated over 
parton distributions. Figure 4 shows the cross section as a function of dijet rapidity difference 
A for fixed and running ag- In both cases the falling parton densities more than compensate 
for the rise in a with the net result that the cross section falls off with increasing A. The 
slower rise in a for the running case that we saw in Fig. 1 translates here into a faster fall-off 
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than one expects in the fixed conpling case. The lowest-order QCD resnlt is shown as a 
solid cnrve. Its slower fall-off than the two BFKL cnrves reflects the fact that it inclndes 
lowest-order kinematics only, whereas the BFKL resnlts inclnde kinematic snppression dne 
to emission of extra glnons. 

Becanse of the snbtleties involved in nntangling the BFKL prediction for a from the 
effects of the falling parton densities in the measnred dijet cross section, it was proposed in 
Refs. P, ^ to measnre the azimnthal angle decorrelation between the two jets. This qnantity 
is relatively insensitive to the details of the parton densities and provides a clear distinction 
between the predictions of next-to-leading-order QCD and BFKL. As noted above, previons 
comparisons were to BFKL predictions which did not acconnt for snbleading effects snch as 
the rnnning of a* and trnncation of transverse momentnm integrals for the emitted glnons. In 
Fignre 5 we present onr prediction for (cos A0) with these effects inclnded. For comparison, 
the analytic BFKL prediction is shown as a solid curve. Clearly, subleading effects are 
quite substantial. In particular we see that the azimuthal decorrelation occurs more slowly 
with increasing dijet rapidity difference than predicted by the analytic BFKL result. This 
can be understood partly in terms of phase space availability — for a given transverse 
momentum threshold, dijets produced at larger rapidities require more energy, leaving less 
phase space available for emission of gluons. This effect - not present in the analytic solution 
- partly mitigates the increasing probability for emitting more gluons in the center-of-mass 
system. Because the additional emitted gluons are responsible for the decorrelation, the 
result of including subleading effects is a reduction in azimuthal decorrelation compared to 
the analytic BFKL solution, as seen in the figure. 

Finally, we return to the question that originally motivated this work: can the inclusion 
of subleading effects improve the agreement between the measured dijet decorrelation and 
that predicted by BFKL? We show for reference in Figure 5 some recent preliminary mea¬ 
surements from the D0 collaboration with the same minimum jet transverse momentum 
of 20 GeV. The comparison should be taken as a rough guide only, because our predictions 
(and our kinematic cuts) are at the parton level, and we assume equal and opposite rapidities 
of the dijets. The D0 measurements (and cuts) are at the jet level, and the net dijet rapidity 
is allowed to range between ±0.5. The finite width of the jets will lead, for example, to 
(cosA0) 7 ^ 1 even at A = 0. Furthermore, the results are not final, and the error bars in 
the figure represent the statistical and uncorrelated systematic uncertainties only; an error 
band showing correlated jet energy scale systematic uncertainties appear at the bottom of 
the figure. See Ref. [0 for a full explanation of the data. Having made those qualifications, 
we note a marked improvement in agreement between data and the BFKL prediction when 
subleading effects are included. 


4 Conclusions 

We have presented the formalism and numerical results for dijet production at hadron collid¬ 
ers in the BFKL approach, using an improved formalism incorporating an iterated solution 
(as described for deep inelastic scattering in [|^) that unfolds the sum over emitted gluons 
that is implicit in the analytic solution to the BFKL equation. We have cast the iterated so¬ 
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lution in the form of an event generator. This allows us to incorporate subleading effects such 
as energy conservation and other kinematic constraints as well as the running of the strong 
coupling constant, which are necessarily absent in the analytic approach. It also allows us 
to examine the properties of dijet events. We hnd that the subleading effects included in 
the improved formalism can be substantial, and in particular they lead to improved agree¬ 
ment with measurements of the azimuthal decorrelation in dijet production at the Fermilab 
Tevatron. Further elaboration of our results for hadron colliders, as well as applications of 
the formalism to forward jet production in ep collisions, will appear in future work. 
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Figure Captions 

[1] (a) The subprocess cross section d for dijet production as a function of the dijet 
rapidity difference A for as fixed (circles) and running (x’s), for minimum jet transverse 
momentum Pt = 20 GeV and /i = 1 GeV. The analytic solution is shown as a solid 
curve. In each case the cross section is normalized to its value at A = 0. Errors are 
from Monte Garlo statistics, (b) The mean of cosA0 computed from the subprocess 
cross section, for the same cases as in (a). 

[2] The subprocess cross section a for dijet production as a function of the gluon resolution 
cutoff /i for A = 1 (circles) and A = 3 (x’s) and Pt = 20 GeV, with running 

[3] The contributions to the subprocess cross section a for dijet production from different 
numbers of emitted (i.e. resolved) gluons ng, for A = 1, 3, 5, running Pt = 20 GeV 
and /i = 1 GeV. 

[4] The total cross section for dijet production at the Tevatron as a function of the dijet 
rapidity difference A for fixed (circles) and running (x’s), for minimum jet transverse 
momentum Pt = 20 GeV and /i = 1 GeV. The lowest order QGD result is shown as 
a solid curve. 

[5] The azimuthal angle decorrelation in dijet production at the Tevatron as a function of 
dijet rapidity difference A for running as (x’s), for minimum jet transverse momentum 
Pt = 20 GeV and fi = 1 GeV. The analytic BFKL solution is shown as a solid curve 
and a preliminary measurement from D0 is shown as diamonds. Error bars on D0 
points represent statistical and uncorrelated systematic errors; correlated jet energy 
scale systematic are shown as an error band. See text for further explanation. 
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